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ABSTRACT 


Topaz  II  is  a  Russian  made  nuclear  powered  satellite  purchased  by  the  United  States.  This 
thesis  provides  a  vibration  analysis  of  the  power  unit  of  Topaz  II.  Several  experimental  vibration 
tests  were  conducted  on  the  actual  power  unit  by  Sandia  National  Laboratories,  in  Albuquerque, 
New  Mexico.  Utilizing  Structural  Dynamics  Research  Corporation’s  I-DEAS  Masters  Series 
Software,  a  finite  element  model  was  made  of  the  power  unit.  Using  the  I-DEAS  software,  this 
model  was  tested  using  the  same  vibration  test  values  as  those  performed  on  the  actual  power 
unit.  The  lateral  and  axial,  experimental  and  theoretical,  sine  tests  were  then  analyzed.  The 
purpose  of  this  analysis  was  to  determine  how  accurate  a  representation  the  model  was  of  the 
actual  Topaz  unit.  This  was  accomplished  through  comparison  of  natural  frequencies  obtained 
experimentally  with  those  predicted  from  I-DEAS  testing  of  the  finite  element  model.  The 
frequency  range  of  interest  was  from  0  Hz  to  200  Hz.  When  solved  for  system  modes,  the  model 
predicted  aU  natural  frequencies  within  4.6  Hz  of  those  found  experimentally.  The  average 
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I.  INTRODUCTION 


A.  ACQUISITION  OF  TOPAZ  U 

The  United  States  conducts  many  of  its  military  missions  in  space  through  the  use  of 
satellites.  The  Department  of  Defense  has  a  long  term  requirement  to  develop  follow-on  power 
technologies  to  support  these  missions.  Having  the  capability  to  produce  a  large  quantity  of 
electrical  power  on  orbit  would  allow  a  greater  spectrum  of  missions  to  be  performed,  such  as 
environmental  monitoring  and  using  high  power  active  sensors  for  surveillance.  To  determine  the 
most  effective  means  of  supplying  this  power,  a  series  of  Air  Force  and  Department  of  Defense 
mission  and  integration  studies  were  begun  in  the  late  1980's.  The  conclusion  drawn  was  that 
compact,  environmentally  safe  nuclear  power  was  the  best  option.  In  the  spring  of  1992,  a 
Department  of  Defense  study  considered  the  option  of  importing  the  Russian  made  TOPAZ  II 
space  nuclear  reactor  technology.  It  was  decided  that  this  technology  would  greatly  enhance 
development  of  U.S.  space  nuclear  reactors.  [Ref  1] 

Procurement  of  TOPAZ  II  was  accomplished  by  the  Strategic  Defense  Initiative 
Organization  (SDIO),  now  known  as  the  Ballistic  Missile  Defense  Organization  (BMDO).  They 
have  also  been  responsible  for  the  transfer  of  space  nuclear  reactor  technology  between  the 
United  States  and  Russia.  This  started  in  December  of  1991,  when  SDIO  began  the  investigation 
of  the  possibility  of  a  U.S.  launched  Russian  TOPAZ  11  nuclear  power  system.  At  this  time,  the 
intended  destination  of  TOPAZ  11  was  in  the  Nuclear  Electric  Propulsion  (NEP)  Space  Test 
Mission.  [Ref  1] 

The  foremost  consideration  in  launching  a  TOPAZ  II  satellite  was  safety.  The 
requirement  was  for  a  safe  launch,  within  constraints  of  the  given  budget.  To  find  out  if  this  were 
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possible,  SDIO  subpoenaed  a  group  of  twelve  engineers  and  scientists  from  Sandia  National 
Laboratories  (SNL),  the  Air  Force  Phillips  Laboratory  (PL),  Los  Alamos  National  Laboratory 
(LANL),  the  University  of  New  Mexico  (UNM),  and  Advanced  Sciences,  Inc.  (ASI).  The 
conclusion  of  this  comprehensive  preliminary  study,  completed  in  1992,  was  that  the  TOPAZ  II 
could  be  launched  safely  and  within  budget  constraints. 

In  the  spring  of  1992,  two  TOPAZ  II  reactors  were  delivered  to  the  New  Mexico 
Engineering  Research  Institute  (NMERI).  The  fuel  elements  were  not  included  in  the  package. 
Also  obtained  by  NMERI  was  a  complete  satellite  qualification  testing  facility.  Under  the 
direction  of  the  Phillips  Laboratory,  non-nuclear  testing  of  TOPAZ  II  was  begun  as  part  of  the 
TOPAZ  Technology  Program. 

The  TOPAZ  Technology  Program,  led  by  the  BMDO,  is  a  three  phase  cooperative  effort 
with  Russia.  The  three  phases  are  comprised  of  joint  evaluation,  joint  testing,  and  joint 
development.  Each  phase  allows  an  opportunity  for  both  the  U.S.  and  Russia  to  evaluate  the 
program  and  determine  future  direction.  The  United  States  stands  to  gain  possession  of 
thermionic  reactor  technology  for  use  by  U.S.  industry.  [Ref  1] 


B.  MISSION  DESCRIPTION 

Mission  goals  of  the  NEP  Space  Test  Program  include  use  of  the  NEP  Space  Test 
Spacecraft  to  test  an  assortment  of  electric  propulsion  thrusters  and  conduct  additional  scientific 
experiments  that  would  be  compatible  with  its  other  mission  requirements.  Two  specific  elements 
of  Russian  space  technology  will  also  be  tested  and  evaluated,  those  being  the  TOPAZ  II  nuclear 
reactor  and  the  SPT-100  electric  propulsion  thruster.  One  of  the  programs  primary  goals  is  to 
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evaluate  the  NEP  in  orbital  transfer  applications  and  measurement  of  the  spacecraft  environment 
produced  by  NEP  system  operation. 

The  reactor  system  of  the  NEP  spacecraft,  shown  in  Figure  1,  is  separated  from  the  main 
section  of  the  spacecraft  by  an  extendable  boom.  This  arrangement  is  necessary  to  maintain  the 
radiation  dosage  received  by  other  spacecraft  components  to  acceptable  levels.  [Ref  1] 

C.  TOPAZ  n  POWER  SYSTEM  DESCRIPTION 

The  power  system  of  TOPAZ  II  consists  of  a  reactor  power  system,  which  generates 
electricity  from  nuclear  heat  through  the  use  of  in-core  thermionic  conversion  units.  This  system 
was  designed  by  Russian  engineers  to  meet  the  following  specifications:  [Ref  1] 

1 .  Under  no  circumstances  should  the  reactor  system  operate  before  reaching  orbit 
altitude. 

2.  The  system  shelf  life  must  be  10  or  more  years 

3.  The  coolant  shall  not  freeze  before  operation 

4.  The  system  must  be  capable  of  providing  6  kWe  at  27  volts  for  a  lifetime  of  3  years 
with  a  rehability  of  0.95. 

5.  The  mass  of  the  power  system  must  not  exceed  1061  kilograms,  excluding  the 
controller  and  supply  system  masses. 

The  main  subsystem  of  the  TOPAZ  II  power  system  are  the  reactor  subsystem,  the  radiation 
shield,  the  primary  coolant  loop,  the  cesium  supply  system,  the  instrumentation  and  control 
system,  gas  systems,  thermal  cover,  and  the  primary  power  system  structure,  (see  Figure  2) 
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Instrument  Boom 


The  nuclear  reactor  is  a  small  zirconium  hydride-moderated,  epi-thermal  reactor  with  in- 
core  thermionic  converters.  The  reactor  incorporates  37  in-core  single  cell  thermionic  fiiel 
elements  (TFEs).  Prior  to  fuel  being  loaded  into  the  TFEs,  electric  heaters  can  be  loaded  into  the 
TFE  cavity  to  simulate  heat  generated  from  the  reactor.  The  advantage  of  using  heaters  is  the 
capability  to  perform  non-nuclear  testing  of  the  thermionic  converters  and  the  complete  power 
system.  With  heaters,  this  can  be  accomplished  at  conditions  near  what  the  nominal  operating 
conditions  would  be. 

The  purpose  of  the  radiation  shield  is  to  attenuate  neutron  and  gamma  radiation.  Located 
behind  the  reactor,  it  is  a  truncated  cone  composed  of  a  stainless  steel  shell  containing  lithium 
hydride  (LiH)  which  is  the  neutron  shield.  The  radiation  shield  is  attached  by  support  legs  to  the 
lower  end  of  the  reactor. 

The  reactor  coolant  loop  is  comprised  of  a  single  EM  pump,  stainless  steel  piping,  sodium- 
potassium  (NaK)  coolant,  and  a  heat  rejection  radiator.  From  a  lower  plenum,  the  NaK  passes 
through  the  reactor  core  where  it  is  heated,  to  an  upper  plenum  to  two  stainless  steel  pipes  to  the 
radiator,  the  NaK  flows  through  two  coolant  pipes  dividing  into  three  pipes  each,  prior  to 
reaching  the  EM  pump.  The  pump  then  sends  the  NaK  back  to  the  reactor  lower  plenum. 

Cesium  is  supplied  to  the  TFE  interelectrode  gap  via  the  cesium  supply  system.  This  is 
necessary  to  increase  the  efficiency  of  the  thermionic  converters  through  suppression  of  the  space 
charge  that  occurs  near  its  emitters. 

The  Instrumentation  and  Control  system  is  required  for  the  monitoring,  controlling,  and 
telemetering  of  the  power  system  status.  Its  primary  duties  are  the  following:  1)  to  start  the 
power  system,  2)  maintenance  of  the  system  operation  under  nominal  operating  conditions,  3) 
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stabilization  of  the  voltage  received  by  the  payload  4)  performance  of  the  commands  from  the 
ground  control  station,  5)  shutdown  of  the  power  system,  6)  maintain  safety  control  during  land- 
based  operations,  7)  telemetering  performance  data  to  the  ground,  8)  shunting  of  excess 
electrical  power  to  ballast  resistors,  and  9)  to  keep  the  storage  battery  charged.  [Ref  1] 

There  are  four  gas  systems  on  board  TOPAZ  II.  These  are  the  COj/He  moderator  cover 
gas,  the  He  TFE-gap  gas,  the  Ar/He  gas  in  the  volume  accumulator  and  the  helium  gas  in  the 
radiation  shield.  The  gas  pressure  is  monitored  during  filling,  electric  ground  testing,  and 
throughout  the  first  few  thousand  hours  of  space  operation  by  pressure  gauges  connected  to  each 
of  these  gas  systems. 

The  purpose  of  the  thermal  cover  is  to  reduce  the  heat  loss  from  the  power  system  on  the 
launch  pad  and  during  the  flight  to  orbit.  If  the  TOPAZ  II  coolant  were  allowed  to  freeze, 
coolant  circulation  would  cease,  and  the  system  would  no  longer  operate  properly. 

The  primary  power  system  structure  consists  of  the  reactor,  the  shield,  and  the  frame. 
These  are  the  three  main  structural  members  of  the  power  system.  The  frame  supports  the 
radiator,  the  gas  bottles,  coolant  piping,  and  the  volume  accumulator.  It  is  attached  to  the  shield, 
and  is  also  used  to  attach  the  spacecraft  boom.  [Ref  1] 

A  more  indepth  analysis  of  the  nuclear  power  system  was  done  by  LT  Steve  Benke  and 
LT  Rich  Venable.  Their  theses  may  be  obtained  through  the  Naval  Postgraduate  School, 
Monterey,  CA.  [Ref  2]. 
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II.  TOPAZ  II  -  THE  MODEL 


A.  DEVELOPMENT  /  PURPOSE  OF  TOPAZ  II  FINITE  ELEMENT  MODEL 
The  Finte  Element  Model  of  TOPAZ  n  was  developed  by  J.M.  Lacy,  of  Idaho  National 
Laboratories  [Ref  3],  One  of  the  first  problems  encountered  was  that  of  which  system  ofxmits 
to  use.  American  engineers  tend  to  favor  the  SAE  system  whereas  the  Russians  use  their  version 
of  the  SI  system  Mr.  Lacy  elected  to  use  the  Russian  modified  metric  system  The  following 
units  conversion  table  may  prove  useful.  [Ref  3] 


Unit 

Dimension 

Equivalencies 

length 

1  mm 

0.03937  in 

0.00329  ft 

mass 

1  kg 

2.2  Ibm 

0.00571  Ib-sVin 

force 

IN 

0.2248  lb 

0.1019  kgf 

stress 

1  MPa 

IN/mm^ 

145.04  psi 

The  finite  element  model  was  made  for  the  reactor  unit  and  its  support  fi-ame.  The  frame 
has  three  legs  whose  feet,  when  connected  by  an  imagmary  circle,  form  a  circle  of 2440  mm  in 
diameter.  The  upper  end  of  the  fi'ame,  to  which  the  reactor  unit  is  attached,  has  a  diameter  of 
680  mm  This  platform  is  2116  mm  above  the  feet,  (see  Figure  3)  On  the  launch  pad,  the  power 
system  and  fi^ame  will  weigh  approximately  1043  kg.  The  center  of  gravity  will  be  2304  mm 
above  the  feet. 

Some  of  the  TOPAZ  fi-ames  are  constructed  of  12X18H10T  stainless  steel 'wMe  others 
are  reported  as  08X18H10T  materiaL  These  designations  are  interpreted  below  and  their 
properties  are  listed  in  Table  1.  [Ref  3] 
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12 

X18 

HIO 

T 

density,  p 

modulus  of  elasticity,  E 


less  than  0.12%  Carbon 
1 8%  Chromium 
10%  Nickel 

Titanium  is  5  times  the  Carbon  content, 
but  less  than  0.8% 

7.95  g/cm'^3 

195.12  GPa  (28.3x10®  psi)* 

*  an  assumption,  to  correspond  with  ASME 
stainless  steel  properties 


Designation 

Temperature 

Cc) 

Ultimate 

Strength 

(MPa) 

Yield 

Strength 

(MPa) 

Elongation 

(%) 

Reduction 
in  Area 
(%) 

08X18H10T 

unknown 

490 

196 

40 

55 

12X18H10T 

350 

549 

216 

35 

NA 

600 

510 

196 

40 

55 

Table  1.  TOPAZ  11  Frame  Properties  [Ref.  3J. 


One  purpose  of  creating  the  FEM,  and  the  one  addressed  by  this  thesis,  is  the 
determination  of  the  natural  frequency  of  the  reactor  system.  This  is  a  major  concern  during  the 
launch  phase,  due  to  the  vibration  forces  that  will  be  imparted  to  TOP.AZ  The  natural  frequency 
determines  how  the  system  is  going  to  respond  to  these  inputs  and  whether  or  not  it  will 
withstand  the  launch  loads.  If  the  natural  frequencies  are  accurately  calculated,  then  the  stiffness, 
mass,  and  mass  distribution  of  the  system  are  known  to  have  been  modeled  correctly.  Having 
accomplished  this,  the  model  can  then  be  used  reliably  to  predict  what  stresses  and  displacements 
will  occur  under  actual  loading  conditions. 

Through  testing  of  an  actual  TOPAZ  unit  in  1992,  the  Russian  scientists  found  the  first 
natural  frequency  of  the  reactor  unit  to  be  4.5  -  5.0  Hz.  This  presented  a  problem  because  the 
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natural  frequency  of  the  entire  system  would  be  far  lower  than  that  of  the  rector  unit  alone.  Due 
to  interference  from  launch  loads,  a  payload  whose  natural  frequency  is  less  than  4  Hz  cannot  be 
launched  in  the  United  States. 

When  E  D.  Schaefer  of  APL  [Ref  8]  created  a  FEM  to  confirm  the  frequency  range  of  4.5 
to  5  Hz,  he  instead  found  the  first  natural  frequency  of  the  system  to  be  closer  to  10.5  Hz.  Early 
in  1993,  Lacy  [Ref  3]  was  tasked  to  study  this  problem  independently.  His  results  were  the  same 
as  Schaefer's.  At  this  point,  the  discrepancy  explanation  was  thought  to  lie  in  differences  between 
the  reactor  system  and  its  documentation.  Through  hands-on  measurements,  this  was  disproved. 
Properties  that  could  not  be  directly  measured  were  confirmed  with  the  Russian  scientists. 

During  the  summer  of  1993,  Sandia  National  Laboratories  performed  a  modal  survey  of 
an  actual  reactor  unit.  This  system  massed  899  kg  (as  opposed  to  1043  kg  used  the  previous 
FEM)  and  the  first  natural  frequency  was  found  to  be  20.3  Hz.  In  their  testing,  Sandia  used  both 
the  reactor  unit  and  the  seismic  mass.  When  this  mass  was  included  in  Lacys'  FEM  and  the 
reactor  mass  reduced  to  899  kg,  he  predicted  the  first  natural  frequency  to  be  21 .6  Hz. 

Therefore,  Sandia's  modal  survey  supported  both  Lacys'  and  Schaefer's  theoretical  results  of 
10  Hz. 

Lacy  produced  another  FEM  of  TOPAZ  in  the  fall  of  1993  [Ref  3].  This  model  weighed 
899  kg  and  had  small  refinements  made  to  the  mass  distribution  of  the  system,  and  to  the  frame 
geometry.  His  results  continued  to  agree  with  Schaefer's  results  and  those  of  Sandia  National 
Laboratories.  The  first  natural  frequency  of  the  seismic  mass  was  only  2%  below  that  reported  by 
the  modal  survey,  and  the  fij'st  lateral  fixed  base  frequency  of  the  system  was  found  to  be  10.8  Hz. 
[Ref  3]. 
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The  results  presented  in  this  work  were  derived  through  vibrational  analysis  of  this 


updated  FEM  model. 
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m.  EXPERIMENTAL  SHOCK  AND  VIBRATION  TESTING  OF  TOPAZ  H 


A.  SETUP 

Early  in  September,  1994,  shock  and  vibration  testing  was  performed  on  TOPAZ  II  at 
Sandia  National  Laboratories  in  Albuquerque,  N.  M.  Specifically,  the  vibration  tests  performed 
were  lateral  and  longitudinal  sine  vibration  tests,  random  vibration  tests,  and  shock  tests.  The  sine 
tests  are  examined  in  this  thesis.  Analysis  of  the  random  vibration  tests  has  been  performed  by 
LT  Sheryl  Campbell  fi'om  the  Naval  Postgraduate  School,  in  a  companion  thesis  completed  in 
June  1995. 

There  were  twelve  accelerometers  mounted  on  the  power  system.  The  accelerometers 
were  placed  where  the  anticipated  responses  were  deemed  most  critical.  Their  locations  are  given 
in  Tables  2  and  3  for  the  lateral  and  longitudinal  tests  respectively.  Also  shown  in  this  table  are 
the  instrument  axes  fi'om  which  responses  were  read  by  the  testing  equipment.  Data  were  not 
taken  from  all  axes  from  each  accelerometer  during  each  test  because  the  testing  equipment  could 
only  read  a  maximum  of  24  channels  concurrently. 

There  were  three  different  coordinate  systems  used  between  the  experimental  set-up  and 
theoretical  analysis.  The  actual  TOPAZ  unit,  the  FEM,  and  the  accelerometers  each  had  their 
own  coordinate  directions.  The  relationship  between  the  instrument  and  accelerometer  axes  is 
given  in  Table  4.  Figure  4  shows  the  relationship  between  the  FEM  and  TOPAZ  coordinate  axes. 
Table  4  does  not  include  accelerometer  numbers  1 1  and  12.  These  accelerometers  were  mounted 
on  the  cesium  unit  and  the  start-up  unit  of  TOPAZ.  In  the  FEM,  these  units  were  not  individually 
modeled  so  that  there  was  no  theoretical  result  obtained  from  these  Jpcations.  In  future  research. 
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Accel. 

No. 


Instrument  Location 


Mount  on  -Z  leg  of  frame  at  base 


Mount  on  leg  between  +Y  and  +Z  of  frame  at  base 


Mount  on  bottom  collector  of  radiator  on  -Z  axis 


Mount  on  bottom  collector  of  radiator  on  +Y  axis 


Mount  on  -Z  leg  of  frame  at  top  of  radiator 


Mount  on  leg  between  +Y  and  +Z  of  frame  at  top  of 
radiator 


Mount  on  reactor  leg  bracket  most  closely  aligned 
with  -Z  axis 


Mount  on  reactor  leg  bracket  most  closely  aligned 
with  +Y  axis 


Mount  on  reactor  top  plenum  on  -Z  axis 


Mount  on  reactor  top  plenum  on  +Y  axis 


Mount  to  "hard  point"  on  cesium  unit 


Mount  to  "start-up  unit"  frame 


Axis  Axis  Axis 
Used  Used  Used 

X  Y  Z 


Table  2.  Lateral  Sine  Sweep  Accelerometers  [Ref.  6]. 


these  units  should  be  modeled  and  their  responses  evaluated. 

In  this  thesis,  the  terms  lateral  and  longitudinal  refer  to  the  lateral  (Z)  and  axial  (X) 
directions  of  the  actual  TOPAZ  unit.  The  position  of  the  accelerometers  on  TOPAZ  are  depicted 
in  Figures  4  and  5.  Also  shown  are  the  coordinate  systems  of  both  the  actual  TOPAZ  unit  and  the 
FEM,  for  future  reference  in  later  chapters. 

The  values  of  the  sine  vibration  applied  to  TOPAZ  were  chosen  as  common  launch 
frequencies  that  might  be  imparted  to  the  satellite.  The  specific  values  and  test  parameters  are 
listed  below.  (A  graphical  repesentation  of  the  sine  input  is  given  in  Figure  6)  [Ref  6] 
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Accel. 

Axis 

Axis 

Axis 

No. 

Instrument  Location 

Used 

Used 

Used 

X 

Y 

Z 

1 

Mount  on  -Z  leg  of  frame  at  base 

X 

X 

X 

2 

Mount  on  leg  between  +Y  and  +Z  of  frame  at  base 

X 

X 

3 

Mount  on  bottom  collector  of  radiator  on  -Z  axis 

X 

4 

Mount  on  bottom  collector  of  radiator  on  +Y  axis 

X 

5 

Mount  on  -Z  leg  of  frame  at  top  of  radiator 

X 

6 

Mount  on  leg  between  +Y  and  +Z  of  frame  at  top  of 

X 

X 

radiator 

7 

Mount  on  reactor  leg  bracket  most  closely  aligned 

X 

with  -Z  axis 

8 

Mount  on  reactor  leg  bracket  most  closely  aligned 

X 

with  +Y  axis 

9 

Mount  on  reactor  top  plenum  on  -Z  axis 

X 

X 

10 

Mount  on  reactor  top  plenum  on  +Y  axis 

X 

X 

11 

Mount  to  "hard  point"  on  cesium  unit 

X 

X 

X 

12 

Mount  to  "start-up  unit"  frame 

X 

X 

X 

Table  3.  Longitudinal  Sine  Sweep  Accelerometer  Locations  [Ref.  6]. 
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SINE- VIBRATION  INPUT 


Frequency(Hz) 


Level:  g  (1  min.  @  axis) 


5 

5-8 

8-40 

40-100 

100-200 


.25 

linear  increase 
1.0 
.9 
.8 


1 .  The  input  was  ramped 

2.  The  test  duration  (sweep  rate)  was  0.25  octaves/minute  which  corresponds  to  a  test 
duration  of  about  7  minutes. 

3.  A  standard  test  tolerances  of  +/-10%  was  placed  on  the  test  parameters. 

4.  Three  additional  control  accelerometers  were  used.  Their  purpose  was  to  limit  the 
response  of  the  TOPAZ  system  to  prevent  overstressing  it.  These  were  accel¬ 
erometers  number  2,  5,  and  9.  See  Figures  4  and  5  for  their  locations.  [Ref  7]. 
These  controls  were  limited  to  responses  of  6,  9,  and  12  dB  respectively.  [Ref  9]. 


The  sine  vibration  was  applied  to  the  TOPAZ  unit  twice,  once  in  the  longitudinal  X 
direction  and  once  in  the  lateral  Z  direction.  Due  to  the  geometry  of  TOPAZ,  the  Z  direction  was 
determined  by  the  Russians  to  be  the  most  critical  of  the  two  lateral  axes  directions  for  testing 
purposes. 
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Accel 

# 

On  Node 

# 

Instrument 

Axis 

1 

5 

X 

1 

5 

Y 

1 

5 

Z 

2 

12 

X 

2 

12 

Y 

2 

12 

Z 

5 

29 

X 

5 

29 

Y 

5 

29 

Z 

6 

55 

X 

6 

55 

Y 

6 

55 

Z 

7 

79 

X 

7 

79 

Y 

7 

79 

Z 

8 

75 

X 

8 

75 

Y 

+Y(30 

degree 

offset) 


+Z(30 

degree 

offset) 


+Y(30 

degree 

offset) 


-X(20 

degree 

slant) 


-Z(20 
degree 
slant  & 
15 

degree 

offset) 


-X(20 

degree 

slant) 


Accelerometer  ID 


I01/C13;-Z  LEG@BASE  Axis:  -X 


I01/C14;.Z  LEG@BASE  Axis;  -Z 


I01/C15:-ZLEG@BASEAxis;  -Y 


I02/C16;LEG/+Y  AND  +Z/BASE  Axis:  -X 


I02/C17:LEG/+Y  AND  +Z/BASE  Axis:  +Y 


102/C  18:LEG/+Y  AND  +Z/BASE  Axis;  -Z 


I05/C3 1 :  JOINT  ON  LEG  Axis:  -X 


I05/C33:JOINT  ON  LEG  Axis:  -Z 


106/Cl 9;LEG  +Y  AND  +Z/RAD.TOP  Axis:  -X 


I06/C20:LEG/+Y  AND  +Z/RAD.TOP  Axis; 
+Z(30) 


I07/C05:REAC.LEGBRACKET/-Z  Axis; 
-X(20) 


I07/C08;REAC.LEGBRAK./-Z  Axis:  -Z(15) 


I08/C0 1  :REAC.LEG  BRACKET/+Y  Axis: 
-X(20) 


I08/C04;REACT.LEGBRAKET/+Y  Axis:  +Z 
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Accel 

# 

On  Node 
# 

Instrument 

Axis 

Topaz 

Axis 

Accelerometer  ID 

8 

75 

Z 

+Y(15 

degree 

offset 

&20 

degree 

slant) 

- - - - - - — - -  - 

9 

97 

X 

-X 

I09/C09;REACT.TOP  PLENUM/-Z  Axis:  -X 

9 

97 

Y 

+Y 

109/Cl  liREACT.TOP  PLENUM/-Z  Axis:  +Y 

9 

97 

Z 

-Z 

109/C  12.  RE  AC.  TOP  PLENUM/-Z  Axis!  -Z 

10 

92 

X 

-X 

I10/C03,REAC.TOP  PLENUM7+Y  Axis:  -X 

10 

92 

Y 

+Z 

10 

92 

Z 

+Y 

I10/C02:REAC.TOP  PLENUM/+Y  Axis;  +Y 

Table  4.  Correlation  of  Coordinate  Systems, 


vi;  o 


Figure  4.  Location  of  Accelerometers,  +XY  Topaz  Plane  [Ref.  7J. 
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Acceleration,  mm/s 


10000 

9000 

8000 

7000 

6000 

5000 

4000 

3000 


B.  DATA  IVIANAGEMENT 


The  data  received  from  these  tests  were  converted  over  to  DOS  readable  ASCII  format. 
Results  were  then  manipulated  by  a  simple  MATLAB  program  and  plotted  as  acceleration  versus 
frequency.  One  such  plot  was  obtained  for  each  accelerometer  axis  from  which  data  were  taken, 
for  both  the  lateral  and  longitudinal  sine  tests.  These  graphs  are  contained  in  Appendix  D, 

Each  of  these  graphs  has  been  analyzed  to  determine  at  which  frequencies  peaks  occurred. 
These  peaks  represent  natural  frequencies  of  the  TOPAZ  system  which  were  excited  by  the 
vibration  input.  Below  100  Hz,  most  peaks  are  easily  distinguishable.  Above  100  Hz,  there  are 
many  smaller  peaks  appearing  in  the  data.  This  was  due  to  a  combination  of  noise  in  the 
equipment  and  interference  from  excitation  of  lower  system  modes.  Only  prominent  peaks  are 
tabulated. 

One  example  of  these  graphs  is  given  in  Figure  7.  Peaks  are  indicated  by  'x'  marks.  The 
corresponding  data  aregiven  in  Table  5.  The  values  were  obtained  by  first  analyzing  the  plot  and 
then  going  into  the  converted  data  files  to  locate  the  actual  frequency  and  amplitude  values.  One 
such  plot  and  table  were  made  for  each  accelerometer  axis  from  which  data  was  taken,  for  each 
test.  The  remaining  graphs  and  tables  may  be  found  in  Appendices  D  and  E,  respectively. 

The  system  was  tested  up  to  200  Hz.  Frequencies  above  this  value  are  insignificant  during 
a  launch  scenario  and  are  typically  determined  during  acoustic  testing.  They  are  also  more 
difficult  to  accurately  distinguish. 
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TRANSDUCER. 


I11/C30:HARD  POINT/CESIUM  UNIT  Axis:  +Z 


PEAK 


1 

2 

3 

4 

5 

6 


FREQUENCY  (Hz) 
7.81651 
24.9373 
60.0392 
71.8386 
97.2413 
135.027 
161.199 


AMPLITUDE  (g) 
2.05755 
9.57841 
0.904517 
0.974768 
0.865449 
0.232790 
8.948593E-02 


Table  5.  Data  Derived  From  Figure  1 , 


IV.  DESIGN  TOOLS 


A.  I-DEAS  SOFTWARE 

I-DEAS  [Ref.  4]  is  a  product  of  the  Structural  Dynamics  Research  Corporation  (SDRC). 
I-DEAS  is  an  integrated  package  of  mechanical  engineering  software  tools.  Among  the 
applications  available  with  the  software  are  solid  modehng  design,  drafting,  finite  element 
simulation,  dynamics,  post-processing  and  testing. 

In  the  vibrational  analysis  of  the  TOPAZ  II  power  system,  the  finite  element  model 
solution,  post-processing,  and  dynamics  applications  are  the  most  useful.  The  model  solution 
application  has  the  capability  of  creating  a  "solution  set"  which  describes  the  type  of  solution 
(statics,  dynamics,  etc),  the  type  of  output  to  store,  and  other  options.  Using  this  application,  the 
TOPAZ  II  finite  element  model  was  analyzed  under  both  linear  statics  and  dynamics.  The 
restraints  used  were  a  clamp  located  on  each  of  the  three  feet,  restraining  the  model  against  all  six 
possible  rigid  body  motions.  The  model  was  solved  in  order  to  obtain  the  natural  frequencies  of 
the  structure.  This  solution  also  gave  the  displacements  and  stresses  for  each  mode.  These 
results  were  listed,  and  displayed  two-dimensionally  and  color-coded,  using  the  post-processing 
appUcation.  The  modes  were  animated  to  determine  each  mode  shape.  Figure  8  is  a  snapshot  of 
mode  one. 

After  the  model  was  solved,  the  dynamics  application  was  used  to  apply  the  vibration  sine 
values  of  acceleration  vs.  frequency  to  the  model.  These  loads  were  applied  through  the  clamped 
feet.  A  viscous  damping  ratio  of  0.02,  determined  from  the  previous  modal  test,  was  applied  and 
the  model  response  obtained  and  graphed.  From  these  graphs  the  theoretical  natural  frequencies 
were  obtained.  [Ref  4] 
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V.  FIMITE  ELEMENT  MODELING  AND  ANALYSIS 
A.  FINITE  ELEMENT  ANALYSIS  -  THEORY 

The  finite  element  method  provides  a  formulation  which  can  use  digital  computer 
automation  for  the  analysis  of  irregular  systems.  To  accomplish  this,  the  method  sees  a  complex 
structure  as  an  assembly  of  finite  elements,  where  each  element  is  part  of  a  continuous  structural 
member.  By  requiring  that  the  displacements  be  computed  and  the  internal  forces  in  balance  at 
points  shared  by  several  elements,  called  nodes,  the  entire  structure  reacts  as  a  single  entity.  The 
finite  element  method  expresses  the  displacement  at  any  point  of  the  continuous  element  in  terms 
of  a  finite  number  of  displacements  at  the  nodal  points  multiplied  by  given  interpolation  functions. 
[Ref  5] 

The  advantage  of  the  finite  element  method  over  other  methods  is  that  the  equations  of 
motion  for  the  system  can  be  derived  by  fiurst  deriving  the  equations  of  motion  for  a  typical  finite 
element  and  then  assembling  the  individual  elements'  equations  of  motion.  The  motion  at  any 
point  inside  the  element  is  obtained  by  means  of  interpolation,  where  the  interpolation  functions 
are  the  same  for  every  element.  Complicated  structures  can  be  modeled  by  simpler  one  or  two 
dimensional  elements,  resulting  in  less  computation  time.  [Ref  6] 

The  I-DEAS  software  utilizes  the  displacement  based  finite  element  method  to  solve  finite 
element  models.  To  derive  the  equilibrium  equations,  the  principal  of  virtual  work  is  used.  The 
principle  of  virtual  work  is  basically  a  statement  of  the  static  equilibrium  of  a  mechanical  system 
and  is  defined  as  follows:  the  work  performed  by  the  applied forces  through  infinitesimal  virtual 
displacements  compatible  with  the  system  constraints  is  zero.  [Ref  5] 
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The  principle  can  be  used  to  calculate  the  portion  of  static  equilibrium  of  a  system.  For  a 


conservative  system,  the  equation  is 


dW  =  F.-6r.  =  ■— 6x.  +  —by.  +  —bz.  =  0 

1=1  i=i  I  ox  dy.  dz.  . 


where 


bW  =  virtual  work 


Fj  =  the  applied  force 


bq  =  virtual  displacement  vector  of  the 
6V  =  virtual  potential 

bxj  =  virtual  displacement  in  the  X  direction 
by;  =  virtual  displacement  in  the  Y  direction 
bzi  =  virtual;  displacement  in  the  Z  direction 


As  utilized  by  I-DEAS,  the  equation  is  as  follows 


l(ef(o)dV  =  (Uf(F.)  +  f(Uf(f^)dS  + 


where 


(  U  )  =  virtual  displacement 


(  o  )  =  strain  due  to  the  virtual  displacement 


( s )  =  actual  stress  at  current  time  t  due  to  applied  loads 
( F,J  =  nodal  forces 
( )  =  surface  forces 
( fg  )  =  body  forces 

This  equation  states  that  the  work  done  by  point  loads,  surface  loads,  and  body  forces  are 
equal  to  the  internal  strain  energy  due  to  the  virtual  displacement  when  integrated  over  the 
volume.  Given  a  set  of  nodal  values,  the  forces  and  displacements  can  be  evaluated  anywhere  on 
the  body.  This  is  accomplished  through  the  introduction  of  (N),  continuous  element  interpolation 
fiinctions,  and  (B),  the  derivatives  of  (N).  Cancelling  the  virtual  displacement  vector  from  above, 
leaves 

Y,  f[Bf(o)dr^  =  (F.)  +  Y  ^  E  (3) 

elan  elan  elan  ^  ^ 

A  global  set  of  equations  are  formed  from  the  integrations  over  each  element.  When  considering 
the  linear  case,  this  equation  is  reduced  to 

[KKu)  =  (F.)  +  (F,)  +  (F,)  (4) 


where 

( u  )  =  displacement 
[  K  ]  =  linear  stifl&iess  matrix 
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The  linear  stiffness  matrix  is  equal  to 


[^]  =  E  {BYmmdV^ 


where 


|  DJ  =  matrix  of  elastic  constants 
To  determine  strain  and  stress,  the  following  equations  are  used 

(e)  =  [5](«) 


(o)  =  [D](€) 


(5) 


(6) 

(7) 


A  finite  element  model  can  be  made  of  more  than  one  type  of  element  The  idea  is  to 
model  the  deflection  of  the  structure,  not  its  appearance.  The  FEM  is  the  idealization  of  the  entire 
structural  problem.  This  includes  the  elements  node  locations,  both  the  physical  and  material 
properties,  boundary  conditions  and  loads.  [Ref  4] 
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1.  Modal  Analysis 

The  FEM  was  tested  using  I-DEAS  software.  The  restraints  used  were  clamps  on  each  of 
the  three  feet.  The  model  was  solved,  using  the  parameters  from  Appendix  A,  to  determine  its 
natural  frequencies.  Because  the  model  used  ideal  clamps,  these  natural  frequencies  were 
expected  to  be  higher  than  those  found  experimentally.  No  losses  were  modeled,  also 
contributing  to  higher  theoretical  natural  frequencies.  These  frequencies  are  listed  in  Table  6. 

When  the  FEM  was  solved,  post  processing  revealed  the  results  shown  in  Table  6.  Using 
I-DEAS,  modes  up  through  200  Hz  were  animated  to  determine  each  mode  shape  present.  The 
upper  bound  of 200  Hz  was  chosen  because  that  was  as  far  as  the  experimental  testing  went.  The 
mode  shapes  demonstrate  how  the  structure  moves  when  a  particular  natural  frequency  is  excited. 
Figures  8  -  12  are  snapshots  of  five  of  these  mode  shapes.  Each  snapshot  gives  the  mode  number, 
corresponding  frequency,  and  color-coded  scale  relating  to  either  stress  or  displacement. 
Displacement  values  are  displayed  in  Figures  8,  9, 10,  and  1 1,  and  stress  values  in  Figure  12. 
Either  parameter  may  be  selected. 

The  numerical  stress  and  displacement  values  on  these  pictures  are  significant  as  relative 
quantities  only.  Using  Figure  9  as  an  example,  it  is  known  that  any  "red"  section  will  deflect  ten 
times  as  far  as  any  "royal  blue"  section  of  TOPAZ.  If  the  figure  were  displaying  relative  stress, 
the  same  would  hold  true  for  values  of  stress. 

When  solving  for  natural  frequencies,  the  I-DEAS  software  is  solving  an  eigenvalue 
problem.  The  shapes  of  the  functions  are  known  as  the  eigenfunctions.  Because  the  equations 
being  solved  are  homogeneous,  the  shape  of  the  function  can  be  detemmined  uniquely  but  the 
ampUtude  of  the  functions  are  arbitrary.  Each  natural  frequency  has  a  corresponding 
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THEORETICAT.  DATA 


MODE 

FREQ 

(Hz) 

DISPLACED 

AMPLITUDE 

RELATIVE 

AMPLITUDE 

ANIMATED 

MODE  SHAPE 

1 

9.2949 

l.OOE+03 

35-36-37 

beam,  structure,  +/-  Y 

2 

9.4301 

1  .OOE+03 

35-36-37 

beam,  structure,  +/-  X 

3 

27.6152 

3.37E+02 

11 

torsion,  about  Z 

4 

40.5769 

8.58E+02 

3 

mode  2  beam,  +/-  Y 

5 

40.9074 

8.46E+02 

4 

mode  2  beam,  +/-  X 

6 

56.7768 

l.OOE+03 

35-36-37 

axial  +/-  Z 

7 

59.0425 

5.33E+02 

6 

bottom  +/-  Y,  top  +/-  Y 

8 

59.05082 

5.36E+02 

5 

bottom  +/-  X 

9 

63.30842 

4.69E+02 

7 

bottom  only,  radial  +/-  Z 

10 

81.85867 

4.29E+02 

9 

structure  mode  3,  XZ  plane 

11 

81.94865 

4.3E+02 

8 

YZ  deflection  of  reactor 
support 

12 

99.2141 

3.55E+02 

10 

torsion,  frame 

13 

115.4492 

2.04E+02 

16 

XY  plane,  deflection  2 
lower  bays,  opposite 

14 

115.6353 

1.71E+02 

21 

mode  3,  ZX  plane 

15 

120.2439 

2.03E+02 

17 

torsion  in  lower  2  bays 

16 

125.5667 

2.99E+02 

13 

axial,  structure 

17 

126.487 

2.44E+02 

15 

mode  3  beam,  XY  plane 

18 

127.1178 

2.69E+02 

14 

mode  3  XZ 

19 

129.8433 

2.02E+02 

18 

torsion  about  Z,  top 
opposite  bottom  half 

20 

130.4791 

1.47E+02 

25 

mode  3  ZY 

21 

130.483 

1.57E+02 

22-23 

mode  3  ZX 

22 

150.6242 

2.00E+02 

19 

mode  3  ZY 

23 

150.7446 

1.90E+02 

20 

mode  3  ZX 

24 

170.473 

3.30E+02 

12 

mode  3  ZX 
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THEORETICAL  DATA.  CONT. 

MODE 

FREQ 

(Hz) 

DISPLACED 

AMPLITUDE 

RELATIVE 

AMPLITUDE 

ANIMATED 

MODE  SHAPE 

25 

174.9176 

1.32E+02 

26 

torsion  of  bottom  bay 

26 

177.8734 

1.28E+02 

27 

mode  3  ZY 

27 

178.1774 

1.48E+02 

24 

mode  4  ZX 

28 

181.3174 

1.24E+02 

29 

mode  4  ZY 

29 

182.0637 

1.15E+02 

31-32 

torsion,  structure 

30 

183.776 

1.57E+02 

22-23 

mode  3  ZX 

31 

185.9247 

8.96E+01 

1 

torsion,  lower  2  bays 

32 

189.7822 

l.OlE+02 

34 

torsion,  structure 

33 

190.003 

8.80E+01 

2 

torsion,  structure 

34 

193.3796 

1.15E+02 

31-32 

axial,  torsion,  structure 

35 

193.6082 

1.26E+02 

28 

mode  3  ZY,  torsion  lower 
half 

36 

199.7701 

1.18E+02 

30 

axial,  torsion 

37 

212.4125 

l.lOE+02 

33 

— 

Table  6.  Theoretical  Results . 
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Figure  10.  Mode  6,  Displacement 


VALUE  OPTION:  ACTUAL 
SHELL  SURFACE:  TOP 
2.  48E  +  06 

2.  23E  +  06 

+ 

1.  98E  +  06 

1.  73E  +  06 

1.  49E  +  06 

1.  24E  +  06 

9.  91E  +  05 

7.  44E  +  05 

4.  96E  +  05 

2.  49E  +  G5 


1-  41E+03 


eigenfunction,  or  natural  mode.  Because  the  problem  is  homogeneous,  any  constant  multiplier  of 
a  natural  mode  represents  the  same  natural  mode.  Therefore,  amplitudes  are  relative  values  only. 
[Ref  6] 

Table  6  also  lists  the  relative  amplitudes  of  the  displaced  mode  shapes.  This  is  important 
to  determine,  when  excited,  which  mode  shapes  will  deflect  the  structure  the  greatest  amount. 
This  information  also  gives  a  feeling  of  where  the  maximum  areas  of  stress  will  occur.  A 
description  of  the  animated  mode  shapes  are  listed  in  the  last  column  of  the  Table.  All  axis 
references  in  this  table  are  to  the  FEM  coordinate  system.  Again,  these  may  be  used  to  predict 
the  structural  deflection  when  the  corresponding  frequencies  are  excited. 
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2,  Forced  Response 

In  order  to  determine  the  accuracy  of  the  FEM,  it  had  to  be  tested  in  the  same  manner  as 
was  the  actual  structure.  Table  7  gives  the  sine-vibration  input  that  was  applied  to  the  actual 
TOPAZ  unit.  This  same  input  was  applied  to  the  FEM.  The  sine  vibration  is  shown  graphically  in 
Figure  6.  This  function  was  applied  to  the  FEM  through  each  of  the  clanged  feet.  A  damping 
ratio  of  0.02  was  applied  to  all  elements  of  the  FEM,  corresponding  to  the  danq)ing  ratio  of  the 
actual  structure. 


SINE  VIBRATION 

Frequency  (  Hz ) 

Acceleration  (  mm/s^) 

5 

2452.5 

5-8 

linear  increase 

8-40 

9810 

40  -  100 

8829 

100  -  200 

7848 

Table  7.  Applied  Sine  Vibration  [Ref.  6]. 


The  sine-vibration  test  was  applied  twice,  once  in  the  lateral  direction  and  once  in 
longitudinally,  as  was  done  in  the  experimental  testing.  The  same  lateral  direction  (Z  of  the  actual 
imit)  was  used  for  testing  of  the  FEM. 

When  post  processing  the  theoretical  results,  data  was  taken  from  locations  on  the  FEM 
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which  corresponded  to  accelerometer  locations  on  the  actual  structure.  As  with  the  experimental 
data,  the  theoretical  results  were  graphed.  Each  of  these  graphs  were  then  compared  to  each  of 
the  experimental  graphs  obtained  from  the  same  location  on  TOPAZ.  An  example  of  one  of  the 
theoretical  graphs  may  be  seen  in  Figure  13,  This  graph  came  from  node  92  of  the  FEM,  which 
corresponds  lO  accelerometer  number  10  of  the  actual  TOPAZ  structure.  This  accelerometer  was 
mounted  on  the  reactor  top  plenum.  (See  Table  3)  As  explained  in  Chapter  II,  this  comparison  is 
necessary  to  determine  the  accuracy  of  the  FEM.  The  results  of  this  analysis  are  explained  in 
Chapter  V. 

In  Figure  13,  under  the  graph,  the  ID  number  of  3021  represents  this  particular  function  in 
I-DEAS.  Each  fiinction  was  assigned  a  different  number  automatically.  The  functions  were  then 
plotted  The  data  from  this  figure  were  obtained  from  the  X  direction  of  the  FEM,  The  Y  axis 
labeled  'modulus'  represents  a  non-dimensional  number.  The  theoretical  amplitude  values  were 
not  necessary  in  determining  which  natural  frequencies  were  excited.  As  with  the  experimental 
data,  prominent  peaks  were  tabulated  from  the  graphs  (See  Figure  13  and  the  corresponding  data 
listed  in  Table  8).  The  phase  plots,  such  as  the  one  in  Figure  13,  were  useful  in  identifying  which 
peaks  were  significant.  When  the  phase  changes  180  degrees,  a  natural  frequency  of  the  system 
has  been  crossed.  As  the  frequency  values  increase,  the  phase  plot  becomes  less  angular  and 
phase  changes  less  distinct.  This  is  because  lower  natural  frequencies  are  also  contributing  to 
phase  changes.  X  marks  are  indicated  on  Figure  13  to  show  which  peaks  were  tabulated. 

The  remainder  of  the  theoretical  graphs  and  tabulated  data  may  be  found  in  Appendices  B 
and  C,  respectively. 
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VI.  ANALYSIS  AND  CONCLUSIONS 
A.  ANALYSIS  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

The  theoretical  and  experimental  graphs  of  Appendices  B  and  D,  respectively,  were 
analyzed  the  same  way.  Prominent  peaks,  indicating  natural  frequencies  of  the  system,  were 
marked  on  the  graphs.  The  raw  data  were  then  examined  so  that  the  frequencies  at  which  these 
peaks  occurred  could  be  exactly  identified.  This  procedure  was  followed  for  both  the  lateral  and 
longitudinal  tests  of  both  the  theoretical  and  experimental  data.  The  lateral  and  longitudinal  test 
results  were  analyzed  separately.  This  was  done  because  different  natural  frequencies  would  be 
excited  during  the  two  tests. 

At  this  point,  there  were  four  sets  of  tabulated  results.  These  were  the  lateral  and 
longitudinal  results  from  the  theoretical  and  experimental  tests.  Each  of  these  sets  were 
individually  studied  to  determine  which  frequencies  had  peaks  commonly  occurring.  The 
experimental  lateral  and  longitudinal  frequencies  are  listed  in  Table  9.  The  theoretical  lateral  and 
longitudinal  results  were  combined  and  listed  as  the  theoretically  predicted  natural  frequencies  of 
the  system.  (See  Table  9) 

As  expected,  the  experimentally  determined  natural  frequencies  were  lower  than  those 
theoretically  predicted.  As  mentioned  earlier,  this  is  in  part  due  to  the  ideal  restraints  used  by  the 
computer  model.  Also,  no  losses  were  simulated  on  the  FEM,  but  losses  would  occur  in  an  actual 
system.  Some  of  the  theoretically  predicted  frequencies  were  not  found  in  the  experimental  data. 
This  was  due,  in  part,  to  the  accelerometer  placements  on  TOPAZ.  They  may  not  have  been 
located  properly  to  register  each  natural  frequency.  Some  frequencies  may  have  been  damped  out 
by  losses  in  the  system.  This  accounts  for  the  "missing"  experimental  frequencies  in  Table  9. 
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THEORETICAL  DATA 


MODE  THEORETICAL  SINE-X  SINE-Z 
NATURAL  ACTUAL  ACTUAL 
FREQ  (Hz)  (Hz)  (Hz) 


9,2949 


9.4301 


27.6152 


40.57694 


44,111 


56.7768 


59.0425 


59,05082 


63.30842 


81.85867 


81.94865 


99.2141 


115,4492 


115.6353 


120.2439 


125.5667 


126.487 


127.1178 


129,8433 


130.4791 


130.483 


150.6242 


ANIMATED  MODE 
SHAPE 


beam,  structure,  +/-  Y 


beam,  structure,  +/-  X 


torsion,  about  Z 


mode  2  beam,  +/-  Y 


mode  2  beam,  +/-  X 


axial  +/-  Z 


bottom  +/-  Y,  top  H-  Y 


bottom  -r/-  X 


bottom  only,  radial  +/-  Z 


structure  mode  3,  XZ  plane 


YZ  deflection  of  reactor 
support 


torsion,  frame 


XY  plane,  deflection  2 
lower  bays,  opposite 


mode  3,  ZX  plane 


torsion  in  lower  2  bays 


axial,  structure 


mode  3  beam,  XY  plane 


mode  3  XZ 


torsion  about  Z,  top 
opposite  bottom  half 


mode  3  ZY 


mode  3  ZX 


mode  3  ZY 


mode  3  ZX 


mode  3  ZX 


150 

169 


154 

170 


THEORETICAL  DATA.  CONT. 


MODE 

THEORETICAL 
NATURAL 
FREQ  (Hz) 

25 

174.9176 

26 

177.8734 

27 

178.1774 

28 

181.3174 

29 

182.0637 

30 

183.776 

31 

185.9247 

32 

189.7822 

33 

190.003 

34 

193.3796 

35 

193.6082 

36 

199.7701 

37 

212.4125 

ACTUAL 

(Hz) 


SINE-Z 

ACTUAL 

(Hz) 


ANIMATED  MODE 
SHAPES 


torsion  of  bottom  bay 


mode  3  ZY 


mode  4  ZX 


mode  4  ZY 


torsion,  structure 


mode  3  ZX 


torsion,  lower  2  bays 


torsion,  structure 


torsion,  structure 


axial,  torsion,  structure 


mode  3  ZY,  torsion  lower 
half 


axial,  torsion 


Table  9.  Theoretical  and  Experimental  Result  Tabulation . 
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The  theoretically  predicted  natural  frequencies  were  then  plotted  against  the 
experimentally  found  frequencies.  This  graph  is  shown  in  Figure  14,  For  70%  of  the  predicted 
modes,  there  exists  a  corresponding  experimentally  determined  natural  frequency.  Modes  where 
experimental  frequencies  were  not  found  were  most  likely  a  result  of  inadequate  accelerometer 
placement.  In  a  one  to  one  comparison,  theoretically  predicted  and  experimentally  found  natural 
frequencies  were  within  4.6  Hz  of  each  other,  with  an  average  difference  of  1 .8  Hz. 


B.  CONCLUSIONS 

The  finite  element  method  is  a  very  useful  tool  for  structural  analysis.  Complex  structures 
can  be  modeled,  have  forces  applied,  and  be  tested  for  failure  in  a  number  of  ways,  using  only 
computer  software.  This  provides  a  substantial  cost  savings  when  used  to  complement  prototype 
or  actual  unit  testing.  The  I-DEAS  software  proved  invaluable  for  the  analysis  work  presented  in 
this  thesis. 

In  order  to  obtain  reliable  useful  data  from  finite  element  modeling,  the  model  must 
accurately  represent  how  the  actual  structure  will  deflect  under  loading.  One  method  of  checking 
the  accuracy  of  the  model  is  to  solve  it  for  its  natural  frequencies.  If  the  system  modes  predicted 
by  the  model  resemble  the  actual  system  modes,  the  model  is  considered  reliable.  Through  the 
testing  completed  for  this  thesis,  the  results  show  the  two  sets  of  natural  frequencies,  theoretical 
and  experimental,  to  be  very  close,  and  therefore  the  model  is  an  accurate  representation  of  the 
TOPAZ  structure. 
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C.  RECOMMENDATIONS  FOR  FUTURE  STUDY 


Further  study  of  the  finite  element  model  test  data  should  be  done  to  locate  the  highest 
amplitude  response,  thereby  identifying  the  frequencies  which  correspond  to  the  worst  case 
loading.  As  mentioned  previously,  individual  components  such  as  the  cesium  block  or  starter  unit 
should  be  analyzed  to  see  if  they  exceeded  the  limits  of  the  test.  If  the  power  system  must  be 
launched  in  an  upright  position  to  conform  to  the  shape  of  U.S.  launch  vehicles,  structural 
analysis  of  the  frame  must  be  done  to  see  if  it  will  withstand  launch  loads.  From  previous  testing 
in  Russia,  the  interface  between  the  shield  and  reactor  can  be  a  problem  source.  Further  research 
needs  to  be  completed  in  this  area. 

When  data  were  analyzed  in  this  thesis,  only  relative  values  of  stress  and  displacement 
were  obtained.  Further  analysis  needs  to  be  done  to  predict  the  actual  stress  and  displacement 
that  will  occur  during  a  typical  launch  sequence. 
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APPENDIX  A, 


This  Appendix  contains  the  parameters  used  in  the  I-DEAS  [Ref.  4]  software  to  solve  the 
finite  element  model. 


Type  of  solution. 

SOLUTION  SET 

Normal  Mode  Dynamics 

BOUNDARY  CONDITION  SET 

Normal  Mode  Dynamics  -  SVl 
Restraint  Set  (Pinned  Feet) 
Constraint  Set  (Base  Constraints) 
Temperature  Set 

Selected 

Selected 

Selected 

Room  Temperature 

OUTPUT  SELECTION 

Mode  Shapes 

Stress 

Store/List 

Store/List 

OPTIONS 

Method 

Mass  Lumping 

Enforced  Motion 

Number  of  Rigid  Body  Modes 

Solution  No  Restart 
Selected 

Selected 

0 

SOLUTION  CONTROL 

Number  of  Flexible  Modes 

Number  of  Iteration  Vectors 
Maximum  Iterations 

Frequency  Accuracy  (digits) 

Start  Vector  Creation 

Frequency  Shift  Method 

Frequency  Shift  in  Hertz 

43 

58 

10 

3 

Automatic 

Shift  and  Number  of  Modes 
0 

A  minimum  of  10  iterations  were  required  to  converge  on  a  solution. 
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APPENDIX  B. 


This  Appendix  contains  54  theoretical  graphs  obtained  using  the  I-DEAS  [Ref.  4] 
software.  The  FEM  was  solved  after  applying  a  sine  vibration  input.  Results  were  taken  from 
nodes  of  interest,  corresponding  to  accelerometer  locations  on  the  actual  unit.  These  results  were 
then  plotted  and  the  graphs  are  contained  herein. 

Pages  54-80  contain  graphs  fi'om  the  longitudinal-sine  input.  Pages  81-107  contain  graphs 
from  the  lateral-sine  input. 
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APPENDIX  C. 


This  Appendix  contains  the  theoretical  results  obtained  from  analysis  of  25  of  the  graphs 
from  Appendix  B.  Only  the  graphs  corresponding  to  experimental  testing  accelerometer  locations 
were  analyzed.  Prominent  peaks  were  located  and  the  corresponding  frequencies  tabulated.  A 
table  was  made  corresponding  to  each  of  the  25  graphs. 

Pages  110-113  contain  results  from  the  longitudinal-sine  input.  Pages  114-117  contain 
results  from  the  lateral-sine  input. 
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APPENDIX  D. 


This  Appendix  contains  experimental  graphs  obtained  from  testing  of  the  actual  TOPAZ 
unit.  Data  were  taken  from  accelerometers  placed  at  various  locations  on  the  unit,  manipualted 
by  a  MATLAB  routine,  and  then  plotted. 

Pages  120-143  contain  graphs  from  the  longitudinal-sine  input.  Pages  144-167  contain 
graphs  from  the  lateral-sine  input. 
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APPENDIX  E. 


This  Appendix  contains  the  experimental  data  obtained  from  analysis  of  the  graphs  of 
Appendix  D.  Prominent  peaks  were  located  and  the  corresponding  frequencies  tabulated.  A  table 
was  made  corresponding  to  each  graph. 

Pages  170-175  contain  results  from  the  longitudinal  sine-vibration  test.  Pages  176-183 
contain  results  from  the  lateral  sine-vibration  test. 
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EXPERIMENTAL  DATA 


TOPAZ  4  X-AXIS  SINE  (LONGITUDINAL) 

TRANSDUCER  A1  :RESP.;INBD.SPIDER  FX.LEG  Axis:+X 


MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

0.979038 

2 

38.5593 

2.55313 

TRANSDUCER  A2:CTL.;TBL.FX.NEAR  INB.SPDR.  Axis:  +X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

0.949486 

2 

40.4287 

1.16645 

TRANSDUCER:  I10/C03:REAC.TOP  PLENUM/+Y  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.25590 

0.949505 

2 

39.1799 

4.51255 
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TRANSDUCER: 


MODE 


I07/C05:REAC.  LEG  BRACKET/-Z  Axis:  -X(20) 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.25590 

0.964197 

40.0062 

6.00587 

111.373 

0.552872 

TRANSDUCER: 

I09/C09:REACT.  TOP  PLENUM/-Z  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.25590 

0.986588 

2 

35.3706 

6.14214 

TRANSDUCER: 

I09/C12:REAC.TOP  PLENUM/-Z  Axis:  -Z 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

39.4811 

0.998568 

2 

51.3925 

0.889658 

TRANSDUCER: 

I01/C13:-Z  LEG  @  BASE  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

0.933774 

2 

39.5012 

2.36357 

3 

118.356 

1.00910 

TRANSDUCER: 

I01/C14:-Z  LEG  @  BASE  Axis:  -Z 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

37.9520 

0.378004 

2 

72.3113 

0.226455 

3 

148.181 

0.398491 
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TRANSDUCER: 


I01/C15:-Z  LEG  @  BASE  Axis:  -Y 


MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

4.42841 3E-02 

2 

15.1196 

0.217778 

3 

39.0735 

1.24001 

TRANSDUCER:  I02/C16:LEG/+Y  AND  +Z/BASE  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

0.956400 

2 

38.0587 

1.68669 

3 

103.69 

0.884210 

TRANSDUCER:  102/Cl 7:LEG/+Y  AND  +Z/BASE  Axis:  +Y 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

0.130565 

2 

39.5012 

0.780018 

3 

64.2438 

0,187713 

4 

83.3868 

0.249129 

5 

154,734 

0.250545 

TRANSDUCER: 

I06/C19:LEG/+Y  AND  +Z/RAD.  TOP  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.26411 

0.947660 

2 

39.5012 

5.00209 

3 

107.438 

0.670597 
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TRANSDUCER: 


I04/C24:BOTTOM  COLL.RAD./+Y  Axis;  +X 


MODE 


TRANSDUCER: 


MODE 


1 


2 


3 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.26411 

0.916813 

38.4494 

6.19911 

112.567 

1.61763 

I03/C27:BOTTOM  COLL.RAD./-Z  Axis;  +X 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.26411 

0.957385 

38.4494 

4.70091 

109.031 

1.38565 

153.015 

0.480957 

TRANSDUCER: 


MODE 


1 


2 


3 


I11/C28;HARD  POEVT/CESIUM  UNIT  Axis;  -X 


FREQUENCY  (Hz) 


7.25590 


39.7964 


123.496 


AMPLITUDE  (g) 


1.08319 


11.0903 


2.02739 


TRANSDUCER: 


MODE 


1 


2 


3 


I11/C29;HARD  POINT/CESIUM  UNIT  Axis;  -Y 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

23.6227 

1.32583 

40.9289 

2.26666 

58.1368 

2.83179 

84.0117 

1.94570 

123.977 

0.913474 

TRANSDUCER; 


MODE 


1 


2 


3 


I11/C30:HARD  POEVT/CESroM  UNIT  Axis:  -Z 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

8.67318 

3.010930E-02 

28.8288 

1.42810 

69.7589 

0.597402 

113.349 

0.901246 

TRANSDUCER; 


MODE 


1 


2 


3 


I05/C31:  JOINT  ON  LEG  Axis:  -X 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.26411 

0.936477 

38.5593 

3.58424 

101.861 

0.693603 

144.553 

0.628134 

■ 

TRANSDUCER; 


MODE 


1 


2 


3 


I12/C34: STARTUP  UNIT  FRAME  Axis:  -X 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.25590 

0.981654 

39.5960 

8.74582 

123.977 

0.874155 

169.133 

0.584549 

TRANSDUCER  I12/C36:STARTUP  UNIT  FRAME  Axis:  -Y(60) 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.25590 

0.114642 

2 

34.0583 

3.08350 

3 

38.9664 

6.18548 

4 

101.493 

0.546350 

5 

172.147 

0.555203 
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Z-AXIS  SINE  (LATERAL)  EXPERIMENTAL  DATA 


TRANSDUCER  ALSEVE  CTL.l;INBOARD  LEG(N.)  Axis:  +Z 

PEAK 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.5904 

0.750334 

2 

28.4307 

0.895611 

3 

169.748 

0.844949 

TRANSDUCER: 

A2:SINE  CTL.2;OUTBOARD  LEG(E.)  Axis:  +Z 

PEAK 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.5904 

0.777016 

2 

28.4307 

0.948742 

TRANSDUCER  I08/C04: REACT.  LEG  BRAKET/+Y  Axis:  +Z 

PEAK 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

24.9373 

1.38918 

2 

81.8337 

0.13952 

3 

85.0912 

0.178205 

4 

89.5111 

0.210022 

5 

96.9556 

r  0.219762 

6 

119.427 

0.206647 

TRANSDUCER  I07/C08  .REACT.  LEG  BRAK./-Z  Axis:  -Z(15) 

PEAK 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

53.5097 

0.440174 

2 

59.7097 

0.318055 

3 

76.1284 

0.176809 

4 

81.2433 

0.269399 

5 

85.0912 

0.214570 

6 

89.8642 

0.268446 

7 

163.492 

5.619752E-02 
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TRANSDUCER: 


MODE 


1 


2 


3 


I09/C09:REACT.  TOP  PLENUM/-Z  Axis:  -X 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

81.4916 

5.806199E-02 

84.7875 

5.008457E-02 

94.4956 

4.543561E-02 

98.3362 

4.003 127E-02 

119.001 

8.006853E-02 

140.005 

8.117978E-02 

158.818 

1.432525E-02 

164.039 

1.186877E-02 

170.054 

7.224624E-03 

171.907 

8.513646E-03 

TRANSDUCER; 


MODE 


1 


2 


3 


109/Cl  1:REACT.  TOP  PLENUM/-Z  Axis:  +Y 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

8.76736 

1.95501 

15.3975 

9.63723  lE-02 

15.7615 

8.896597E-02 

28.3918 

2.16219 

54.1565 

0.603319 

59.8494 

0.212064 

69.5546 

0.141099 

84.2027 

0.188112 

111.574 

2.786758E-02 

118.594 

2.285344E-0Z- 

143.565 

2.197792E-02 

190.176 

1.028353E-02 
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TRANSDUCER:  I09/C12:REACT.  TOP  PLENUM/-Z  Axis:  -Z 


MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

9.40651 

4.61358 

2 

29.4789 

3.97788 

3 

53.8231 

1.36627 

4 

100.092 

0.138697 

5 

133.096 

5.688959E-02 

6 

175.661 

1.350727E-02 

TRANSDUCER: 

I01/C13:-Z  LEG  @  BASE  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.59040 

0.548643 

2 

24.5878 

9.492084E-02 

3 

28.0407 

0.114665 

4 

32.1489 

8.86981 8E-02 

5 

46.3899 

4.471 105E-02 

6 

56.3866 

5.088181E-02 

7 

68.8677 

4.105160E-02 

8 

81.9021 

1 .659305E-02 

9 

85.1832 

2.000362E-02 

10 

120.925 

0. 140088 

11 

138.978 

0.265800 

12 

177.79 

0.110414 

TRANSDUCER: 

I01/C14:-Z  LEG  @  BASE  Axis:  -Z 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.76132 

0.270319 

2 

10.5609 

1.10978 

3 

143.486 

1.13443 
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TRANSDUCER: 


MODE 


1 


2 


3 


I01/C15:-Z  LEG  BASE  Axis:  -Y 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

6.31161 

1.959258E-02 

8.56065 

6,319125E-02 

12.5204 

5.286855E-02 

28.0407 

0.291 127 

138.978 

0.518210 

TRANSDUCER:  I02/C17:LEG/+Y  AND +Z/BASE  Axis: +Y 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.71202 

0.502636 

2 

28.53 

0.625009 

3 

60.5565 

0.728032 

133,629 


0.683575 


TRANSDUCER: 


MODE 


1 


2 


3 


I02/C18:LEG/+Y  AND  +Z/BASE  Axis:  -Z 


FREQUENCY  (Hz) 


7.59040 


26.8651 


131.239 


AMPLITUDE  (g) 


0.807592 


1.05284 


1.16256 


TRANSDUCER: 


MODE 


1 


2 


3 


I06/C19:LEG/+Y  AND  +Z/RAD.  TOP  Axis:  -X 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

8.28011 

0.625170 

28,4307 

0.209106 

44.6078 

0.180537 

59.7707 

0.154713 

81.3292 

0.109272 

100.051 

0.152750 

161.436 

8.821841E-02 
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TRANSDUCER  I06/C20:LEG/  +Y  AND  +Z/RAD,  TOP  Axis;  +Z(30) 


MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.75082 

1.56896 

2 

24.1758 

2.22697 

3 

28.8881 

2.56046 

4 

55.9942 

0.349077 

5 

96.9556 

0.466277 

6 

136.480 

0.899289 

TRANSDUCER:  I04/C23:BOTT.  COLL.  RAD./  +Y  Axis;  +Y 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

9.86712 

0.105801 

2 

27.943 

2.34252 

3 

48.6237 

4.5289 

4 

69.3226 

2.20957 

5 

88.3314 

3.07370 

6 

130.777 

2.95019 

7 

154.658 

2.88126 

TRANSDUCER  I03/C25;BOTT.  COLL.  RAD./  -Z  Axis;  -Z 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.71202 

0.866554 

2 

27.2555 

1.15438 

3 

48.2127 

1.22134 

180 


TRANSDUCER  11 1/C28:HARD  POINT/CESIUM  UNIT  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

8.76736 

0.281597 

2 

13.1312 

0.125737 

3 

28.8881 

5.36258 

4 

59.9416 

0.819486 

5 

81.4916 

0.468706 

6 

140.465 

0.537631 

7 

192.935 

0.201526 

TRANSDUCER:  I11/C29:HARD  POINT/CESIUM  UNIT  Axis:  -Y 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

8.83778 

1.23716 

2 

25.1265 

5.35297 

3 

38.5431 

1.97241 

4 

60.1536 

3.74441 

5 

87.5979 

0.767632 

6 

97.6267 

0.604192 

7 

135.027 

0.270747 

TRANSDUCER:  I11/C30:HARD  POINT/CESIUM  UNIT  Axis:  +Z 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.81651 

2.05755 

2 

24.9373 

9.57841 

3 

60.0392 

0.904517 

4 

71.8386 

0.974768 

5 

97.2413 

0.865449 

6 

135.027 

0.232790 

7 

161.199 

8.948593E-02 

6 

7 


TRANSDUCER; 


MODE 


1 


2 


3 


I05/C31:  JOINT  ON  LEG  Axis:  -X 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.76132 

0.410090 

28.53 

0.241601 

51.9497 

0.122066 

73.0626 

0.256050 

144.013 

0.365256 

TRANSDUCER: 

I05/C33:JOINT  ON  LEG  Axis:  -Z 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.71202 

1.01215 

2 

27.2555 

1.40077 

3 

140.974 

2.00604 

TRANSDUCER;  I12/C34: STARTUP  UNIT  FRAME  Axis:  -X 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

8.26743 

0.459574 

1 

2 

29.4789 

0.831606 

3 

53.8231 

0.136480 

4 

106.063 

0.593274 

TRANSDUCER: 


MODE 


1 


2 


3 


I12/C35:  STARTUP  UNIT  FRAME  Axis:  -Z(60) 


FREQUENCY  (Hz) 

AMPLITUDE  (g) 

7.75082 

1.83180 

27.8961 

4.99722 

55.9942 

0.397595 

132.123 

0.834946 
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TRANSDUCER:  11 2/C36: STARTUP  UNIT  FRAME  Axis:  -Y(60) 

MODE 

FREQUENCY  (Hz) 

AMPLITUDE  (g) 

1 

7.97612 

1.13312 

2 

24.2569 

1.43728 

3 

28.3918 

1.85856 

4 

59.4898 

0.371952 

5 

108.367 

0.362733 

6 


139.493 


0.220601 
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